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Using the general solution to the Einstein equations on intersecting null surfaces developed by Hayward, we investigate the non- 
linear instability of the Cauchy horizon inside a realistic black hole. Making a minimal assumption about the free gravitational 
data allows us to solve the field equations along a null surface crossing the Cauchy Horizon. As in the spherical case, the 
results indicate that a diverging influx of gravitational energy, in concert with an outflux across the CH, is responsible for the 
singularity. The spacetime is asymptotically Petrov type N, the same algebraic type as a gravitational shock wave. Implications 
for the continuation of spacetime through the singularity are briefly discussed. 
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I. INTRODUCTION 

What are the generic features of gravitational collapse? 
We do not, as yet, have a complete answer to this ques- 
tion, however most relativists are confident that gravita- 
tional collapse leads to the formation of black holes (at 
least within the astrophysical context). Furthermore it is 
plausible that the external field of such a black hole set- 
tles down to a member of the Kerr- Newman family, since 
these are the unique stationary solutions of the electrovac 
Einstein equations (see for example |^). Externally, devi- 
ations from these solutions are expected to die away with 
an inverse power-law in advanced time leaving only the 
mass, charge and angular momentum observable. A ques- 
tion which one might hope to answer is whether such a 
property continues to hold inside the black hole - that is, 
does the internal geometry also approach Kerr-Newman 
form? It seems not 

The known exact black-hole solutions possess Cauchy 
horizons - null hypersurfaces which are the boundary 
of the future domain of dependence for Cauchy data 
of the collapse problem. These horizons exhibit highly 
pathological behaviour; small time-dependent perturba- 
tions originating outside the black hole undergo an in- 
finite gravitational blueshift as they evolve towards the 
Cauchy horizon (CH). This blueshift of infalling radia- 
tion gave the first indications that these solutions, which 
so well describe the exterior geometry at late times, may 
not describe the generic internal structure. Penrose ||] 
pointed this out more than twenty five years ago, and 
since then linear perturbations have been analysed in de- 
tail 1^. These observations led to the conjecture that a 
scalar curvature singularity would form either at or before 
the CH once back-reaction was accounted for. 

Poisson and Israel ||] , generalizing work of Hiscock , 
have shown that a scalar curvature singularity forms 
along the CH of a charged, spherical black hole in a 
simplified model. This singularity is characterised by 
the exponential divergence of the mass function with ad- 
vanced time. The key ingredient producing this tremen- 



dous growth of curvature is the blueshifted radiation flux 
along the CH, although it is also necessary that some 
transverse energy flux be present. In Q it was argued 
that the physics underlying the analysis was sufficiently 
general that similar results should hold for generic col- 
lapse (i.e. upon relaxation of the assumption of spher- 
ical symmetry). Further calculations support the con- 
jecture that the singularity inside a generic black hole is 
nuU §0. 

The aim of this paper is to present a detailed analy- 
sis of the CH inside a non-spherical black hole by taking 
advantage of the recent result due to Hayward He 
showed how to obtain the general solution of the Einstein 
equations on a pair of intersecting null surfaces. We begin 
by recasting the results of Poisson and Israel in a sug- 
gestive form which more closely resembles the approach 
taken to the general problem. Our purpose in section ^ 
is to stress the main points which must be taken as as- 
sumptions in the later calculations. We also feel that 
the analysis highlights the non-linear nature of the effect 
and the merely precipitous role played by the outfiux in 
this model . With this preliminary review of the mass- 
inflation phenomenon out of the way, section HI begins 



the analysis of the CH singularity inside a realistic black 
hole. The assumptions which we make about the na- 
ture of the CH are discussed, and following Hayward |Tl| ] 
we formally integrate the first order Einstein equations 
of appendix A. Using these results we can show that 
the divergences which arise are integrable, and obtain 
the leading behaviour of all quantities of interest near to 
the CH. The asymptotic expressions for the Weyl scalars 
show that the spacetime is Petrov-type N near the CH, 
the same algebraic type as a gravitational shock wave. 
These results (which are also suggested by the previous 
work once again raise questions about the classical 
continuation of spacetime through the singularity. 

The equations and curvatures are relegated to the ap- 
pendices in an effort to maintain the clarity of the pre- 
sentation. Appendix A gives a summary of the notation 
and lists the dual null Einstein field equations in their 
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first order form, as derived by Hayward . Appendix B 
lists the components of the Riemann tensor necessary to 
analyse the algebraic type of the spacetime. 



II. SPHERICAL BLACK HOLE INTERIORS 

In this section we review the mass-inflation phe- 
nomenon in the context of charged, spherical black holes 
as studied by Poisson and Israel The presenta- 

tion differs slightly from that in the literature and 
closely parallels the method used later to discuss the non- 
spherical problem. In this way the limitations of the later 
analysis, and of the approximations used should be more 
apparent. 

The physics behind mass-inflation is relatively simple. 
The CH inside a Reissner-Nordstrom black hole is a null 
hypersurface corresponding to infinite external advanced 
time. Time-dependent perturbations which originate in 
the external universe are gravitationally blueshifted as 
they propagate inwards near to the CH. Thus a charged 
black hole which deviates only slightly from Reissner- 
Nordstrom in the exterior is expected to have a barrier 
of radiation, with an exponentially diverging energy den- 
sity, streaming along parallel to the CH. Generically some 
of this ingoing radiation scatters off the gravitational po- 
tential inside the black hole, leading to a flux of energy 
crossing the CH . It is the non-linear gravitational in- 
teraction of these two fluxes of radiation which generates 
a divergence of the local mass function. It is important to 
realise that the gravitational blueshift of time-dependent 
perturbations is the key ingredient producing this mass- 
inflation. 



The Poisson-Israel model 

It is convenient to use null coordinates on the "radial" 
two spaces so that the spherical line element is 



2e ^dudv + r''{dx'' + sin^ xdy^) 



(2.1) 



where A = X{u,v), r = r{u,v) and the coordinates are 
such that u is a retarded time and v an advanced time. 
The stress-energy tensor for a radial electromagnetic field 
is 



Ef," ^q'^/Snr^ diag(-l, -1, 1, 1) 



(2.2) 



where q is the electric charge on the black hole. Pois- 
son and Israel used crossflowing null dust to model the 
perturbations of the geometry, arguing that the large 
blueshift near to the Cauchy horizon should make the 
Isaacson effective stress energy desciption valid for 
the ingoing radiation. They also pointed out that the 
nature of the outflux is not important, its only purpose 
is to initiate the contraction of the Cauchy horizon. The 
stress-energy tensor is therefore 



where = — and n 



Poutnf_i,ni, , (2.3) 
= —d,,u are radial null vec- 



tors pointing inwards and outwards respectively, and, pi„ 
and Pout represent the energy densities of the inward and 
outward fluxes. Each term in Eq. ( ^.3[ ) is independently 
conserved so that 



L\n{v) 



Pout 



(2.4) 



The functions iin(u) and iout('^) are determined by the 
boundary conditions, however, it is important to note 
that they have no direct operational meaning since they 
depend on the parametrization of the null coordinates. 

The field equations can now be written in a first order 
form by defining the extrinsic fields 



(2.5) 
(2.6) 
(2.7) 
(2.8) 



where a comma denotes partial differentiation. Along S 
there are two evolution equations 



'-2v),u = —{r'-M') 



(2.9) 
(2.10) 



and a focussing equation which describes the behaviour 
of the null cones with vertices at r = 0: 



(r^ e)^u + r^i^- 0/2) 9 = -2Lout(u) 



(2.11) 



The complete set of Einstein equations, including those 
which hold on 5"+ , may be obtained from the above equa- 
tions (p.9[), (2.10) and (2.11) via the symmetry operation 



(u; 9, v, 0, v) {v; 0, 9, v) 



(2.12) 
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It is convenient to imagine that the inflow (outflow) is 
turned on at some advanced (retarded) time v = vq {u = 
uo)- In the pure inflow regime the spacetime is described 
by a charged Vaidya solution 

ds^ = dw{2dr - fdw) + r^{dx^ + sin^ xdy^), (2.13) 



where / is given by 



/=1- 



2m(w) 



(2.14) 



and w is the standard external advanced time coordi- 
nate. In particular it is infinite both on future null in- 
finity and the CH. The only non-zero component of the 
stress-energy tensor is 



2 



T — 



1 dm 



(2.15) 



In the outflow region the solution is of the same form 
with advanced replaced by retarded time. The structure 
of the spacetime is shown in Fig. 0. 




FIG. 1. A spacetime diagram showing the spherical Pois- 
son-Israel model. The influx (outflux) of lightlike dust is 
turned on at advanced (retarded) time vq (uq). The surface 
S~ coincides with the Cauchy horizon, CH. S'^ is parallel to 
the event horizon, EH, and intersects the Cauchy horizon at 
S, on which / = and r — To- 



The solution 

As mentioned above our treatment differs from that in 
the literature We consider the solution to the Ein- 

stein equations only on the two intersecting null surfaces 
S~ and . By choosing S~ to coincide with the CH 
we can show that a scalar curvature singularity must be 
present on the CH in the cross-flow region. 

We first integrate the equations on the CH. This is 
readily achieved by making an appropriate choice of the 
retarded time co ordi na te u , in (2.1^ such that v = 0/2 
along S'~. Then (2.8), (2.11) and (2.6) can be integrated, 
in that order, to give 



49" -2 diiLontiu) I 



rt^ri^riB" u-2 / du \ du Lout{u) 
Jo Jo 



(2.16) 
(2.17) 

(2.18) 



The (sub-) super-script (— ) is used to indicate that these 
equalities hold only on S~ , while an (o) indicates the 
value of the function on the two-sphere S ^ O . 



Of the remaining two equations we only need (|2j 
gives 



du ( - i 



which 



(2.19) 



Eq. ( 2.10D could be integrated in the same way. 

The same construction works on 5"*" (we will use this 
as our primary tool in the non-spherical case) however it 
is more convenient to work with the exact solution ( 2.13| ) 
to determine 9°, 9^ and tq here. 



Assumptions 

The two essential features of this model will be adopted 
in the non-spherical case with only slight modification. 
We wish to emphasise them at this point. 

The first assumption is that there should exist a sta- 
tionary portion of the CH in the spacetime. This is 
achieved by turning on the outflux at some retarded time 
(uq) inside the event horizon of the black hole (see Fig. 
|l|), thus set 



t(u) = (3H{u ~ Uq) 



(2.20) 



where /? is a constant, and H{u) is a Heaviside step func- 
tion. This is probably the most contentious issue in the 
Poisson-Israel ^ and subsequent analyses . It presup- 
poses that the CH begins at finite radius, and essentially 
assumes that the singularity at the CH will be null. Some 
arguments have been advanced to suggest that this issue 
is important p^ , however no convincing evidence has 
emerged to refute its correctness. Indeed preliminary nu- 
merical results [15 1 support this assumption, in contrast 
to the claim in |1£]. 

The second assumption pertains to the influx of radi- 
ation; it is taken to decay with an inverse power law of 
advanced time. This was initially based on an extrapola- 
tion of results due to Price ^ to the event horizon of the 
black hole. Recently it has been verified numerically 
This decay is correctly reproduced by the ansatz 



m(w) 



(2.21) 



where Kq is the surface gravity of the inner horizon, a 
is a dimensionless constant which depends on the lumi- 
nosity of the collapsed star, and p > 12 is an integer 
(the numerical value derives from Price's analysis which 
shows that the radiative tail of the collapse decays as 
an inverse power p > 4Z + 4 of advanced time, where 
I is the multipole order of the perturbation). Both lin- 
ear and non-linear |^ perturbation analyses suggest 
that this inverse power-law decay of perturbations is also 
correct near the CH inside the black hole. 
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Mass inflation 
Along 5*+ the radius obeys the first order equation 



dr 
dw 



1 



1 



2m{w) 



(2.22) 



where m{w) is given by ( p^ . We are interested in cal- 
culating 



= 2r"^ — — 
^ dw dv 



(2.23) 



where v is the coordinate which appears in (2.1) chosen 



so that A = In r along 5+ . It is a straightforward m atte r 
t o ca lculate dw/dv provided one notes that Eqs. ( |2.5|) - 
(2.11) are form invariant under functional rescalings of 
the coordinates u and v. Therefore r also satisfies 



d"^ 2 ^ ( In 
dw"^ dw 



dw 
dv 



d , dm 



on 5+ in view of ( pl2| ) and ( pll| ). Solving ( |2^ ) for r, 
(2.24) can be integrated to get 



^(1 + ...) 



(2.25) 



as w — !■ oo. Thus 



r+ ^ r„ + '^P^' { 1 + (p - 1)(- In \v\)-^ + . . .} 



(p- 1)Ko 



(2.26) 



and 



2a(-ln|w|)"P 



{l+p(-ln|«|)-i + ...} (2.27) 



as w 0. This result means that, while the radius of 
the two-sphere S is finite, 6° is actually infinite (in these 
coordinates). 

Using the invariant definition of the mass in a spherical 
geometry 



2m(x") 



(2.28) 



we know that 6'+6'+ ^ as w ^ 0, which means that 
e° = 0. Substituting ( ^.20|) into ( ^.17|) and using ^° = 



61" = -2(m - uo)(3H{u - uq) . 
The radius of the CH is then 



(2.29) 



c- = a 



(2.31) 



from which we can immediately see the reason for the 
Cauchy horizon singularity. The product 9^6^ — at S 
(see Fig. ^ and hence C~ is finite, however, if /3 is non- 
zero C" jumps to infinity at u — uq when the outflux 
is turned on [see Eq. (2.9)]. Furthermore Eq. (2.3C) 
shows that the radius remains finite for some amount of 
retarded time, so it must be the mass function which 
becomes infinite at u = uq. 

This analysis emphasises that the nature of the out- 
fiux is largely irrelevant, it simply serves to initiate the 
contraction of the CH. It is this contraction which is at 
the root of mass-inflation |q| . 



III. NON-SPHERICAL BLACK HOLES 

Charged spherical black holes have served as a useful 
model in which to investigate the non-linear instability of 
the CH inside a black hole, however, the question needs 
to be addressed in the more realistic non-spherical con- 
text. Here one expects that a black hole settles down 
to a member of the Kerr-Newman family at late times. 
(Our analysis focuses on such a situation when there is 
no electromagnetic field.) 

The Kerr solution can be written as 

y\ sin'^ 9 
ds^ = T;{dw-ad(j)sm^ 6)"^ H ^{V^ + a^]d(j)~adwY 



+2dr{dw-ad(l)sin'^ 0) + p^dO^ 



(3.1) 



where A = — 2 



mr 



9 9 9 



a 



^ cos^ 9 and 



w 



is standard advanced time. The CH, located at r = 
m — \Jw? — and w = oo, is a stationary null hyper- 
surface - its lightlike generators have zero rate of expan- 
sion. Moreover, these generators are shear free since Ray- 
chaudhuri's equation implies that shear produces con- 
traction. This is not a generic situation; a linear per- 
turbation analysis shows that there is usually a flux of 
backscattered radiation crossingthe CH, along with the 
blueshifted influx parallel to it Q. Based on this obser- 
vation and using the spherical case as a guide, we will 
show that the CH is the locus of a scalar curvature sin- 
gularity. The leading divergence may be associated with 
propagating modes of the gravitational field, manifesting 
itself in the form of a diverging Weyl curvature of Petrov 
type-N. 



= rl - {u- uo)^(3H{u ~ Uq) , 



(2.30) 



and finally 9 is given by (2.19). The square of the Weyl 
curvature on the CH is 



Formalism 

We employ the dual null formalism of Hayward |l^] 
which is based on a 2-1-2 null decomposition of space- 
time iflSl. In terms of coordinates u and v which label 
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the intersecting null hypersurfaces that foliate the space- 
time, the line element can be written as 

ds^ = -2e^^dudv + {saS°')du'^ + 2sadx''du + habdx^'dx'' 

(3.2) 

where the lower-case latin indices range over {1,2}. 
There remains freedom to rescale u and v, and to trans- 
form the coordinates so that s" vanishes on a chosen 

V = constant surface, we take this to be the CH located at 

V = (see Fig. ^). hab is the 2-metric, s° a shift 2-vector 
and A a scaling function on S{u, v), the spatial 2-surface 
which is the intersection of the null surfaces labelled by u 
and V. The 2-metric can be decomposed into a conformal 
factor f2 and a conformal 2-metric kab so that 



fl ^hab 



(3.3) 



Thus kab has unit determinant. Each of the quantities 
defined here depend on all four variables u, v, a;". 

The vacuum Einstein equations ( A7)-( |A23 ) are written 
in a first order form in terms of the fields 



^n-^Cin, 

(Tab — ^Clkab, 
V = ClX, 

1 



'^ab = ^ ^nkab 



- -e^VLkabCis^ , 
where C indicates the Lie derivative, and 



(3.4) 
(3.5) 
(3.6) 

(3.7) 



I - — n - — - s" — 

dv du dx°- 



are null vectors. The notation and the equations are 
given in appendix A, more detail may be found in Hay- 
ward |l|]. 

Our aim is to calculate the Weyl curvature scalars on 
the intersecting null hypersurfaces 5'+ and S~ by making 
an ansatz for the free gravitational data {kab) on these 
surfaces. S~ is taken to coincide with the CH at t; = 0, 
while S*"*" is an outgoing null hypersurface, parallel to the 
event horizon, crossing S~ near to P in Fig. |^. The choice 
of the data is motivated by our understanding of the 
spherical situation and by the non-linear perturbation 
analysis of Ori S] . 



Assumptions 



The focussing equation on 5"+ is [Eq. (AT)] 



(3,9) 



Hayward noticed that making use of coordinate freedom 
on 5+ (rescaling v), this equation may be linearised. 
Thus choosing v such that 



(3.10) 



it is possible to obtain the solution to the Einstein equa- 
tions on 5*+. Each of the equations ( A7)-( A15|) can be 
integrated (in the order they appear) provided one knows 
kab on this hypersurface, since 



= (fcf 9.fc+)(fcf ) 



(3.11) 



We cannot evolve generic initial data from the event 
horizon of a non-spherical, rotating black hole so we make 
an ansatz for this conformal 2-metric. It is inferred from 
the work of Ori [|| that the conformal metric may be 
written as 



k+ 



kUx^)+K+,[\n\vlx-) 



where det(fc°^) = 1 and asymptotically 



K 



ah ' 



Tab{x'){- \n\v\y 



(3.12) 



(3.13) 



as w — > 0. The functions J-ab{x'^) are well behaved for 
all a;", and n is an integer. (The value implied by Ori's 
work is n > 6.) This form is based on the observation 
that non-linear metric perturbations decay according to 
an inverse power-law of advanced time |^ , and the expec- 
tation that Hayward's coordinate is related to external 
advanced time by 



(3.14) 



as w — > (or w oo). While this cannot be proven to 
hold, it seems reasonable as long as the surface does 
not encounter a caustic {fl — 0) at or before the CH - 
this condition is satisfied a postiori implying, at least, a 
self-consistent treatment. 

As in the spherical analysis, we give no detail of the 
gravitational data on 5*^ save to point out that the grav- 
itational energy flux crossing the CH should not diverge. 
In fact, on physical grounds, it is expected that any radi- 
ation from the star will be exponentially redshifted near 
P in Fig. ^. Furthemore once the incoming radiation gets 
scattered by the gravitational potential inside the black 
hole, it should only produce a slow contraction of the 
CH ]l| 

To summarise, we assume that a portion of the CH ex- 
ists which is caustic free near to P in Fig. ||. Furthermore, 
the gravitational perturbations propagating into the hole 
decay according to an inv erse p ower law of advanced time 
which is related to v by ( 3.14 ). 
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FIG. 2. A schematic representation of the generic black 
hole interior. The null surface S~ coincides with the Cauchy 
horizon, CH. The other null surface S"*" is at it = and inter- 
sects the Cauchy horizon at 5 before a caustic forms. Initial 
values of the fields are specified on the intersection between 
the hypersurface v = vo and 5*^ . Also indicated is the ingoing 
gravitational wave tail, and its scattered component 



Solution on S 

With this aspect in hand we can proceed to the solu- 



tions of Eqs. (A7)-(A14). Integrating (3.9) 



(3.15) 



A (sub) super-script o indicates the initial value of the 
function at in (see Fig. ||). On this two surface we 
expect all the fields to be well behaved, bounded (and in 
some cases non-zero) functio ns of . 

Based on the assumption ( 3.12 ) and ( 3.13| ) about fc^^ 
we can estimate the behaviour of this integral near to the 
CH. Firstly note that the inverse of the conformal metric 
can be written 



l.ab _ i,ab 



(3.16) 



where e is the 2-dimensional, antisymmetric matrix 
with e^^ — 1. The functional form of K'^j^ implies that 



1 dK+ 
u i91n \v\ 



will diverge since 



(9 In \v\ 



"-?^ab(- In I w| 



-n-l 



(3.17) 



(3.18) 



as u — > 0. Since (Tab<^°'^ is w ^ times a slowly varying 
function (near the CH), 9^ is asymptotically 



ft ^hc^ da 

4w(-ln|i;|)2("+i) 



(3.19) 



which diverges at the CH for the same reason that dyK^f^ 
does. 

It is now straightforward to formally integrate the re- 
maining equations, however our main objective is to show 
which quantities diverge and which remain bounded at 
the CH. Observing that the integral 



i-^{-\n\i\)-'''di 



1 



(-ln|i;|)-"+^ (3.20) 



n ~ 1 

is finite as z; ^ (provided n> 1), it is clear that 

1 



n — 1 



J^ah(-ln|'y|) 



-n+l 



(3.21) 



as ti — > 0, and remains bounded at the CH. Each of the 
equations must be treated in turn, examining the diver- 
gent part of the source and its integral near the CH. We 
quickly see that fl is bounded at the CH, and so all the 
intrinsic quantities on the two dimensional surfaces of fo- 
liation remain finite all the way to the CH. To see this, 
first integrate (A8) to get 



A+ = Ao - ^ / dv 



(3.22) 



then substituting in the asymptotic form for wc find 
'dve+ ~ :^kfkfTi,cT,a fdv ^-i(-ln|«|)-2("+i) 



4(2n-H 1) 



and h ence A+ constant as w ^ 0. The gauge choice 
( |3.10 ) implies that 



(3.24) 



and so il and kab are both bounded at the CH. Now 
continuing the integration 

















Jvq 








[ dv 







A, 



and 



(3.25) 

"AfcA+) , 

(3.26) 
(3.27) 



dv 



r, „-3A4 



+ -i(2)^ + c^>+-iA"A,A+ 



-iA'^A+A„A+ +cj;A,A+ - ^a^l J , 
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from which it is clear that only finite terms appear in the The integration now proceeds as it did for S~^, and in 
equation for 6, which therefore cannot diverge on CH. particular we arrive at expressions for d~ and cr^j, 



a+ ~ (finite) - ^6^+/^+ / dv e-^+eh'^^aab , (3.28) 



which is finite despite the presence of divergent terms 
in the integrand [whi ch lo ok like v~^x inverse powers 
of (— Inji)!), see Eq. ( 3.20|) ]. This analysis implies that 
6*+ and cr^j, diverge inside the black hole, and the diver- 
gence is of the same nature as in the spherical model once 
<Jab<^°'^ is interpreted as the gravitational energy crossing 
S+. 



IV. THE CAUCHY HORIZON SINGULARITY 

The solution obtained in the previous section can be 
used to show that a scalar curvature singularity is present 
on the CH due to the blueshifted influx of radiation. On 
S~ the shift vector is set to zero by a coordinate trans- 
formationi" = 5;°(m, x^) without altering the form of the 
metric ( |3.2| ). (These coordinates exist provided s° i s reg - 
ular on S~ in the original coordinate system. Eq. ( 3.2(: ) 
suggests that this is true.) It is assumed that this trans- 
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this in no way alters 



formation was made in section 
the previous analysis. 

Hayward's scheme to integrate the equations on S' 
can now be completed. Making the gauge choice 



(4.1) 



Eq. ( Al6| ) imphes that 



i;=0 



du {kl'duk^Jik'L^duk^J . (4.2) 



cd ' 



As mentioned in the introduction Q (and all other quan- 
tities) are formally known once fc~j(w, a;"^) is specified. 
Our analysis assumes that k~^ is a slowly varying function 
of u which tends to a non-degenerate value as u ^ — cxd; 
that is, near to P (in Fig. |2[) kah should approach the 
value obtained on the CH in Kerr spacetime. Further- 
more 0~ — > and the spacetime should be asymptotically 
Kerr in the limit m — > — oo. It is not a generic situation 
to have 0~ = and/or = on the CH, since some 
radiation (both from the star, and from backscattering 
of the gravitational wave tail inside the black hole [| 12|) 
will usually be present. This discussion implies that there 
should exist a caustic free portion of the CH near to P 
on which our analysis is valid. 



•0=0 



(4.3) 



the shear of the ingoing null rays orthogonal to S{u = 
0,v), is also bounded at the CH. Finally 

^+=^o+ fdi {ia+af - \e+e+ (3.29) ^nd 



„2A_ 



du 



r.„-3X/ 1(2) 



'R 



'db 



(4.4) 



-ih^a^^' / du {e-^6h'Vadb~ ^e~''^h'V {u:(dUJb) 
Jo 

- 5A(rfAf,)A + iA(rfAAt)A - uj(d^b)^ + ^{d^^b)}) 

- Kb^'- I 

Jo 

- tj^AcA + Acw'^) . 
We saw in the previous section how 



du e"^^i 



UJ UJc 



iAcA'=A + iAcAA'A 



lim 



9+ 



oo 



(4.5) 



thus 9~ and cr"^ are generically unbounded on the CH. Of 
the remain ing fi elds, only v contains divergent quantities. 
From Eq. (A23) we see that 



" 1 1 

du her .a"'' - -6 



(4.6) 



which is actually infinite on CH. Imposing the field equa- 
tions on the Weyl scalars one finds that all but ^'4 contain 
divergent quantities, and hence the CH is singular. 

It is possible to say a little more about the nature 
of the singularity by examining the rate of growth of 
the curvature alo ng S '^. Once again, imposing the field 
equations in Eq. ( Bll ) gives 



1 „2A^a^ 



{2A0'ab + 2Wab - 0'am'l™"0'„b}|_ 



Using cr^jj, etc. the asymptotic expression is 



(finite) X {v^ 



Inlvl 



\n+l 



y 



(4.7) 



(4.8) 



as w — > 0. This is the contribution from Lia^^^, which is 
damped by the smallest number of powers of (— In |ti|). In 
a similar manner the leading behaviour can be extracted 
from ^'i and '^2 giving 

4-+ ~ (finite) X {z;(-ln|t;|)"+i}^^ (4.9) 
~ (finite) X {v{-ln\v\)"+^y'^ (4.10) 

where these leading contributions arise due to the pres- 
ence of terms involving the shear (cr|^^) of the surface . 



7 



A little more work shows that is finite on the CH; no- 
tice that eq. ( ^.26| ) combined with the choice s" ^ on 
S~ , and knowledge that A"*" and w+ are finite on 5*+ n 
implies that s+ ~ w as w ^ 0, hence s!pE+„s" — > on 
5^. Finally ^'4 is independent of 6*+ and cr^^ and there- 
fore is finite. These results imply that the square of the 
Weyl tensor is dominated by ^"0, 



~ (finite) X {w2(-ln|i;|)"+i}"^ + 



(4.11) 



as V 



0. 



V. DISCUSSION 

The treatment of the non-spherical gravitational col- 
lapse is generally very difficult, however it seems likely 
that the external field settles down to that of a Kerr 
black hole. Inside the black hole some progress may also 
be made by advancing two assumptions: 1) There exists 
a caustic free portion of the CH, which is a null hyper- 
surface. 2) There is an ingoing gravitational wave tail 
which gets blueshifted by an infinite amount at the CH 
(this is achieved here by our ansatz for fc^^). Indeed us- 
ing Hayward's formalism the nature of the CH can be 
investigated in some detail. 

we 



Following from the analysis in sections III and fV 



conclude that the CH is the locus of a scalar curvature 
singularity, and the curvature is an integrable function of 
the advanced time. Ori has argued that the latter 
property of the curvature suggests an extension of the 
spacetime through the singularity may be possible. It is 
worth commenting a little further on this point. We have 
demonstrated the asymptotic values of each of the Weyl 
scalars, in particular it has been shown that contains 
the leading divergences. The algebraic classification of 
the Weyl curvature is obtained by solving the quartic 



(5.1) 



for a and exa min in g the degeneracy of its roots 



Based on eqs. (4.8)-( 1.10 ) there is a four-fold degeneracy 
as u — > implying that the gravitational field is asymp- 
totically type N, with repeated principal null direction 
k = n + am + a*m -I- aa*l (where a ^ as v ^ 0). 
It is worth noting that gravitational shock waves are of 
this algebraic type (and are the only curvatures which 
may be confined to a thin skin without a surface layer 
of matter being present |^^); therefore the intuitive pic- 
ture of the singularity which emerges is that of a gravi- 
tational shock propagating along the CH. These results 
might be taken to support arguments for the existence of 
a classical continuation of the geometry beyond the CH 
singularity. However, this is a highly speculative point 
and it should be noticed that the Weyl curvature con- 
tains contributions which are not present in an exactly 
type N spacetime. Indeed these contributions diverge at 



the CH a nd are manifest in the square of the Weyl ten- 
sor ( 4.11 ) (which vanishes for a pure gravitational shock). 
Thus a classical continuation is unlikely, as pointed out 
in , since non-classical matter is needed to confine the 
divergent curvature to a thin layer along the CH. 

To conclude, the hairy singularity inside a generic black 
hole is dominated by the propagating gravitational wave 
tail of the collapse (which decays in the exterior of the 
black hole). 
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APPENDIX A: NOTATION 



In this appendix we summarise the dual nuU formahsm of Hayward pi[ . The reader is referred to his articles for 
more technical details. This approach, based on a 2+2 decomposition, assumes that the spacetime can be foliated 
locally by compact orientable 2-surfaces S. Therefore given a Lorentzian manifold M (with signature — \~ ++) we 
assume that there exists a smooth embedding : S* x [0, C/) x [0, V) M with the induced metric, hab, on S 
being spatial. Latin indices range over {1,2}. In ter ms of a basis which is Lie-propagated along d/du and d/dv the 



spacetime line element may be written as in Eq. (3.2). The Lie derivative along a vector ^ is denoted The natural 
covariant derivative of the metric hab is A^, and the Ricci scalar is denoted by '^^^i?. Tensor valued quantities on the 
2-surface S are denoted X°' ' b... (sometimes we will include a superscript (2) to emphasise the 2-dimensional covariant 
nature). Introducing the null vectors 

Z''^ (1,0,0,0), n^ = (0,l,-s'^). (Al) 

we define 

Sab = C,lhab , Sab = Cnhab ■ (A2) 

To write the Einstein field equations in a first order form Hayward introduces the extrinsic fields: 

6 = i/l°''I]af,, 9 = ^h'^^'Sab, 

C^afc = Safe — Ohab, ab = ^ab — Shab, (A4) 

V ^ Ci\, V = CnX , (A5) 

= ^e'^KbCis'' . (A6) 

These quantities are readily interpreted geometrically; [9, 9) are the expansions in the null directions normal to the 
foliation 5*. The tensors {aab,Oab) are the shears, which are traceless by definition, the vector uja is the twist. The 
two functions {v^v) are called the inaffinities and measure the non-affine quality of the coordinates u and v. Then 
the Einstein equations and contracted Bianchi identities are: the I equations: 

Cl9 =-i02_^^_ l^^^^ah 

Ci\ = V (A8) 

LiVL ^9n (A9) 

Cikab = n-^<yab (AlO) 

ClLUa ^-9uJa + ^A'aba + ^^a{iy-9)-^9Aa\ (All) 

Cis" = 2e-^Lo"- (A12) 

£i9 =-6i0 + e-^(-i(2)i^ + ^^^'^_ iA"AaA + iA„AA"A + tj'^AaA-Aaw'') (A13) 

^lO'ab = '^ach'^'^ l^db + \d (^ab ~ (Tab 

+ 2e'^{ujaUJb - iA(<jAb)A + iA(<jAAb)A + w(aAb)A - A(<jt^b)) 

- e-^ hab {i^cLo" - 5 A'^AcA + i AcA A^A + uj^A^X - A^w^) (A14) 

CiT^ - iffabCT"'- ^^e^+e-^-i(2)i? + 3w^c^'^_ lAaAA-'A-w^AaA) (A15) 

and the n equations: 

£ j = _ :;0^ i^^^^ab ^A16) 

£„A = V (A17) 

jCri^ = 9Vt (A18) 

Cnkab = n-^'^ab (A19) 

CniJa = -euJa - ^A^'^ba - ^Aai^ ^9) + ^9Aa\ (A20) 

£nd =-99 + e-^{-^^^^R + UJaUj''^^A''AaX+^AaXA''X-Uj''AaX + AaUj'') (A21) 
^n(^ab ~ <^ac h"^"^ db + hS<^ab ~ h&^^ab 
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+ 2e ^{uJa(^b - ^A(aAfc)A + iA(aAAfc)A - W(aAb)A + A(aW,,)) 



hab {i^cuj" - iA'AcA + i A,A A"A - w'A^A + A, 



CnV 



jCTabO- 



-A 



(-i(')i?^ + SLJaOj'' - iA,A A"A + Uj'^AaX) 



(A22) 
(A23) 



APPENDIX B: RIEMANN TENSOR AND WEYL SCALARS 



1. Riemann Tensor 



Here we list all the components of the Riemann tensor for the metric (|3.2|): 



^vuv 



avb 



R. 



avb 



-Cl V - s''AaV + ^e^S^Y.abS^V + 2s''Ll UJa - e-^uJaA^X + 'Se-^uj^uja + s^Y.abto" 
+ ie^s"^ {Cl ^ab) + i.s"Eafc A^A - ie-^ A'^A A,A - ie^s'^I],fc/^''^Eeds'^ 



L,i Zjab I- 2'- VZjab ^ 46 ^am'^ ■^nb 

1 „A 



^aub 
^vab 

R-avb = 5e'^£„ Sab + 5e'^S™(AmEafc — A(,Sam) + 5 AaAfeA — -J AqA AftA — AbLOa 



abd 



2^ ^ao"^ "^c I 4 
1 _A „m V 1 



1„A 

4 



1 ^rn y 



CnCJa + ys"' (S"A„)E, 



1 



e^E„,/i™"E„,s« + i A,!? + i (s"A„) A,A 



is™A„, A A„A + ie^s"S„„, s"A„A 



_1 „m 
2 



1 „A 

4 



+ 26 (S LUn) S h,na~(S LUm) UJa + 2^ ^7na + S A^UJa 

^e^^abv + ie^s" A,„A Eat + ie^£„ E^fc + ie^s™ A^E^b - A^A s^E^ 

1 o-^ ('c ^,,™^ V Ic-^V umnyy I 1„A / 1 



= i/i="(AaE,„fc - A„E„fc) - h-''sJ^^R"\an - le\s^Cn ^ab ' ie^s^s^A^E^fe - is^AfcA.A 

+ \e^s'AbX S"'^a^ + ^e^s's"'Ab^am + s'AbLUa + ^S^A^A AfcA - ^S^LObAaX + lAaX h"''^„,b 



ie^s™Ea„,s^«;, - -e^s"Ea„/i^"E„fc - ^s^cJaAbX + s^uOaUJb - \h'"^^„,b^a 



1 „m A \ V _ I 1 ,-,A 



1 ATcm ny y 



1 ^-'^ „" 



e^E„,ft(i;i™s"E„,„ + e-^w;^ 



2^ 



A^A) 



(Bl) 
(B2) 
(B3) 
(B4) 
(B5) 
(B6) 



(B7) 



(B8) 
(B9) 



e s"E(j„/i'^™E;,m 



(BIO) 



2. The Weyl Scalars 



We choose a complex-null tetrad {e'^Z, n, m, m} such that 2m^'^m'^'> — h^^'^ and, Z and n are given by (3^). In this 
tetrad the five Newman-Penrose Weyl scalars are in a Ricci flat spacetime 

*o = je^\2Ci J^ab + 2;/E,b - E„„/i""Eb„)m'^m^ (Bll) 
^-1 = le\-2AaV + 2c^"E„„ + 4A t^a + A^AE^Jm'' (B12) 

^-2 = -j{2e^Cn ^ab + 2AaAbX - AaXAbX - AAbLUa - e^E„,b/l™"Ea„ + 2UaAbX 

+ 2ujbAaX - AuaLOb + e^AfcAs"Ea„ + 2p/s"E„„cc;b)m'^m^ (B13) 
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a 



(B14) 
(B15) 
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